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Abstract 

The Holstein-Primakoff mapping for pairs of bosons is extended in order to accommodate single 
boson mapping. The proposed extension allows a variety of applications and especially puts the 
formalism at finite temperature on firm grounds. The new mapping is applied to the 0{N + 1) an- 
harmonic oscillator with global symmetry broken down to 0{N). It is explicitly demonstrated that 
./V-Goldstone modes appear. This result generalizes the Holstein-Primakoff mapping for interacting 
bosons as developed in ref. M. 
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1 Introduction 



Boson expansion theory (BET) has played a significant role over the past decades in our understanding 
of the nuclear many-body problem. Starting with the pioneering work of Marumori and co-workers 0] , 
and of Belyaev and Zelevinsky the interest in this subject has culminated in the eighties through 
the formulation of the interacting boson model Q. Of particular interest, not only for the many-body 
problem, but also, as it has become apparent, for quantum-field theory (QFT), is the perturbative 
boson expansion (PBE) approach. Extensive use of it has been made in nuclear physics, in order to 
extract anharmonicities beyond the Random-Phase Approximation (RPA) (see ref.|5|, |(| for reviews). 
Up until very recently its application to QFT has not attracted much attention and, therefore, has not 
been fully developed so far. The Holstein-Primakoff mapping for boson pairs, first introduced in pi H], 
was recently applied, however, to the O(N) vector model JlJ. It was demonstrated that the mapping is 
able to systematically classify the dynamics according to the 1/iV-expansion, rendering a promising and 
efficient alternative to the well-known functional methods. Furthermore, considering the model in the 
phase of spontaneously broken symmetry, the powerful machinery of the PBE approach as developed 
for deformed nuclear systems could be transcribed to QFT. As a consequence the Goldstone theorem 
as well as the whole hierarchy of Ward identities were exactly satisfied |jj . 

However, the PBE in general, and the Holstein-Primakoff mapping (HPM) in particular ||, rely on 
the bosonisation of pairs of particles. Thereby, images for particle pairs are generated in an ideal 
Fock-space, while single particle images are absent after mapping. This problem has been appreciated 
for the fermionic case in the early days of the boson expansion theory. Marshalek has proposed an 
extension of HPM for fermions in order to allow for a perturbative boson expansion for both even and 
odd nuclei [1C]. In the present letter, we point out the occurrence of the same problem in the case of 
the PBE for purely bosonic models. 

The need for an extended bosonic HPM to include single bosons clearly revealed itself in fij where the 
lack of ideal single-boson states was an obstacle for defining unambiguously the two-point function for 
the Goldstone mode. While to leading order in the 1/A-expansion this problem was circumvented in 
|lj , a next-to-leading order calculation makes an extended version of the HPM mandatory. 
Finite-temperature applications of the PBE approach is another issue where an extended HPM to 
include single bosons is definitely called for. In the following, we wish to sketch a derivation of an 
extended version of HPM for bosons. We will also discuss an application to the 0{N + 1) anharmonic 
oscillator where it will be demonstrated explicitly that this new method is capable of including single- 



boson images with the correct asymptotic energy. 



2 Extended Holstein-Primakoff Mapping 

As a starting point let us consider a system with two types of bosonic creation and annihilation 
operators: a + , a, and 6 + , 6. Pairing these in all possible ways leads to ten group generators of the 
non-compact Sp(4) group. The pairs a + a, aa, a + a + and analogously the pairs of 6-operators form 
two commuting Sp(2) subgroups. The number conserving bilinears a + a, 6 + 6, a + 6 and ab + span a 
closed U(2) algebra. There remain the bilinears a + b + and ab which do not belong to any non-trivial 
subgroup of Sp(4). 

Our goal will be to first set up the boson images of the ten group generators, replacing in the end the 
6-operators by c- numbers (the condensate). This will lead us to the boson image of the semidirect 
product group Sp(2) (g) -/V(l) made up of the elements a + a, aa, a + a + , a, a + , and 1^ respectively. The 
latter is the desired system because it involves even and odd number of boson operators. 
We will follow earlier work for interacting Fermions by Evans and Kraus Klein, Cohen, Li, 



Rafelski, and Rafelski [12, 13| in which a mapping for the ten generators of the SO(5) group was 
derived. Especially, use is made of the work of the latter group of authors to derive this time the 
mapping of the ten generators of the non-compact Sp{4) group mentioned above. Since there is no 
room to go into details (which will be presented elsewhere) we essentially will only give the result here. 
One first realizes that the six generators of the two commuting Sp(2) algebras can be mapped via the 
usual HPM. The difficult task lies in finding an adequate mapping for the generators a + b + and ab, 
which allows one to close of the full Sp(A) algebra. The reader is invited to consult reference [^] for 
a similar derivation. Introducing a set of three new bosonic operators a, A\, and A 2 , one can show 
that the net result for the complete mapping reads 



a + a + )i = AfJ2 + 4(ni+m), (aa)j = ((a + a + )i) , (a + a) I = 2n 1 +m 



(6+6+)/ = A+J2 + 4(n 2 + m) , (66)/ = ((6+6+) 7 ) + , {b + b) I = 2n 2 + m 



[a + b + )i = a + y2 + 4(m + m)^J2 + 4(n 2 + m) $(m) + A^(rn)AjA{a , 
{ab), = ((a+6+) / ) + , 
(a+b)! = ~ [(bb) I ,(a + b+) I ] , (6+a), = ((a+6),) + , (1) 

where ni,rt2, and m are occupation number operators defined by 

n i = A+A i (i = l,2). (2) 



The +-sign in the Holstein-Primakoff square-root indicates the non-compact character of the group 
at hand. Finally, the function $ is given by 

_ i 

,2 



3>(m) 



r + m 



(3) 



4(m + l)(2m + l)(2m - 1 
where r is a constant which is fixed using physical conditions as will be discussed in the next section. 
These results constitute only an intermediate step towards our final goal. As stated in the introduction, 
one wishes to extend the usual HPM for boson pairs, in such a way as to allow the mapping of single 
bosons, as well. In other words, and following the original Belyaev-Zelevinsky approach, one needs to 
achieve a realization of the following algebra 



[aa , a + a + ] 
[aa , a + a] 
[a , a + a + ] 
[a , a + a] 



2 + 4 a+a 
2 aa , 
2a+ , 



(4) 



where all other possible commutators are assumed but not explicitly shown here. This is nothing but 
the algebra of the semidirect product group Sp{2) ® N{1). The first two commutation relations in 
Eq. (B) remind us of the Sp(2) algebra, and as such, one can propose the bosonic HPM as a second 
realization for it. Here again, the difficulty lies in finding an adequate mapping for the single bosons 
so as to close the algebra above. A way out is to notice that, by considering the limit in which the 
operator b and 6+ are transformed into the identity operator, one can ultimately contract the whole 
Sp{A) group to a non-isomorphic semidirect group Sp{2) ® N(l). This singular transformation, which 



can be thought of as a contraction a la Inonu-Wigner or Saletan [14, 15], gives a clear hint on how to 
proceed with the desired extension. Indeed, the single bosons can be deduced from the contraction of 
the generators a+6+ , ab, a + b, and ab + down to the generators a and a+. With this intuitive picture in 
mind, one can show that the following mapping for the five relevant generators constitutes a realization 
of the algebra in Eq. (|j) . 

(aa)j = \J2 + 4(ni + m) A\ , (a + a + )j = (aa)f , (a + a)i = 2n\ + m , 

(a)/ = ^2 + 4(m +m)Ti(m)a + 2a + AxTx(m) , (a + )/ = (a) + , (5) 

where the occupation number operators are, as before, given by n\ = AfA± , m = a + a, while the 
function T\ reads 



Ti(m) 



z\ + m 



2(m + l)(2m + l)(2m- 1) 



(6) 



Here too Z\ is a constant which will be fixed by using physical conditions as will be explained in the 
next section. It is straightforward to verify, through a direct evaluation of the commutators in Eq. (f|), 
that this is indeed a proper realization. This completes our considerations concerning the mapping. 
In the next section, the formalism will be applied to the interesting case of ./V oscillators and used to 
develop the 1/iV-expansion. 



3 The 0(N + 1) Anharmonic Oscillator 



As an application, let us consider the anharmonic oscillator with an 0(N + 1) symmetry broken down 
to 0(N). The properly scaled Hamiltonian of the system is given by 
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(7) 



Here we have considered an explicit [r\ ^ 0) and a spontaneous ((X a ) ^ 0) symmetry breaking along 
the X G mode. The variables X n ,X a and their conjugate momenta P^^P a are expressed in second 
quantization as 



X-n 



x n 



1 



2u> 
1 



(a + a 4 ") , 
+ b + ) , 
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Pa 



(8) 



The frequency, £ a , of the mode X a will be fixed later. The subscripts ir and a are used in analogy 
with the linear cr-model in QFT, where these modes represent the pion- and sigma fields respectively. 
To sort-out the dynamics according to the 1/A^-expansion, one needs to adapt the mapping derived 
in the previous section to the situation of N oscillators. This can be done in a straightforward way. 
It can be shown that the mapping in this case takes the form 



(aa) i 
{ai)i 



2N + 4(ni +m)Ai , 



(a + a)i = 2n\ + m , 



(a + a + )l = (aa)} , 
(afh = (a t )f , (9) 



'2N + 4(m + m)T N (m)ai + 2 ctj Ax T N (m) , 

where ./V is an integer, ri\ = AfA\, and m = Y^i a t a i> while is a generalization of the T± function, 
of the last section to the case of N oscillators. It reads 

Zn + m 2 + m(N — 1) 



T N (m) 



(10) 



2(m + l)(2m + N)(2m + N - 2) 
The constant will be fixed below. One can also easily verify that this mapping leads to a realization 
of the following algebra: 



[(aa) , (a+a+)l 



2N + 4 (a+a) 



[(aa) , (ct^a)] = 2 (aa) . 
[a i ,(a+a+)] = 2a+, 
[a* , (a+a)] = a* . (11) 

For a finite N, the 0(iV + 1) anharmonic oscillator is purely quantum mechanical. For an infinite 
number of degrees of freedom N — ► oo, on the other hand, it can be used to mimic the quantum- field 
situation of the breaking and restoration of a continuous symmetry. 

Using the mapping in Eq. (^), one can expand the Hamiltonian of the system in powers of the 
operators A, a, b and their hermitian conjugates. One then arrives at an expansion of the form 
H = H (0) + + W (2) + H (3) + W (4) + where the superscripts indicate powers of operators 
without normal ordering. This expansion is in fact not unique and therefore the preservation of the 
symmetries is not necessarily guaranteed. A more useful approach consists in organizing the expansion 
in powers of the parameter N, such that H = NHq + y/NHi + Hi + -7= H3 + jjH^ + ... 
This is possible if one chooses a coherent state as the variational ground state for the model 



exp 



(A 1 )A+ + (b)b + |0>. (12) 



This trial vacuum state must accommodate two condensates respectively for the X a mode and the 
newly introduced boson A\ (see ref. [jl|] for details). The mode a, on the other hand, is not allowed 



to condense. The ground-state energy, NHq = ^J,?)^ , calculated on the coherent state, takes the 



following form 

Ho = ^ (2^ + l) + 9 — (d+ VTTa^) 2 + (d+ VYTd^Y + ^ + gs' - V s , (13) 



where we have introduced for convenience the rescaled condensates s = -7= (X a ) ,d = J j* (A). 



The coherent ground state is fully determined by requiring that the values taken by the two condensates 
above lead to the minimum of Hq. The minimization procedure with respect to s and d gives the 
following two coupled BCS equations 



+ 4gs 2 + — (d+ Vl + d 2 



LO V / S 



2ujdVl + d 2 + ( d + v 7 ! + d 2 ) A = 0, (14) 



where A = _|_ Jj^ _|_ ^/l + d 2 ^j is the gap parameter . 

To gather the full dynamics of the leading order in the 1/iV-expansion one needs to generate the terms 
Hi and H2 of the Hamiltonian. This can be done by using the parameter differentiation techniques 



(see ref. [|J for details). The net result for both Hi and H2 then reads: 
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(15) 



Here, Wo is a constant. Since one is not particularly interested in the ground-state energy, the latter 
will not be further specified. The shifted operators A\ = A\ — (A±), (3 = b — (b), and h\ = A^A\ 
annihilate the coherent state The frequency £ a of the X a mode is fixed such that the bilinear 
part of H2 in the (3 operators is diagonal. It is purely of perturbative character, and the frequency is 
explicitly given by 



£l = i o 2 + l2gs 2 + ^-(d+Vl + d^) . 
Using the gap equations (|l4|) and the following easily verifiable identities 



A = -2dVl + d 2 



CO 



d- VTTd 2 



uj + A = (1 + 2d 2 ) 



(d- Vl + d 2 



(16) 



(17) 



one can establish that Hi vanishes at the minimum. From H2, and more precisely from the coefficient 
of m = J2fLi Oifai, one can deduce the existence of ./V uncoupled modes. The common frequency of 
these modes is denoted by S n and is given by 



£l = to 2 (d - Vl + d 2 ) 4 = cj 2 + igs 2 + 



2g 
f 



V 
s 



(18) 



These N modes are our first asymptotic states. Furthermore, it can be easily verified that they have 
Goldstone character. In other words, their frequency vanishes in the exact symmetry limit (77 = 0), 
and for a finite condensate (s ^ 0). 

It is evident from the ansatz above that the model suffers from infrared divergences. However, since it 
is used for demonstration purposes only, we choose to disregard this difficulty here. Clearly the new 
and important result that has been obtained shows up in the fact that the proposed mapping provides 
asymptotic states in the ideal Fock-space which correspond to the images of the single bosons. It 
should be stressed that this is a non trivial-finding which, as shown above, is a direct consequence of 
the extended HPM. It reproduces the result anticipated in the introduction, in clear departure from 
the HPM for boson pairs |j] and in accordance with the Goldstone theorem. 



So far, only the mapping of the bilinears in Eq. (|9|) was involved in expanding the Hamiltonian. The 
single-boson part of the mapping, on the other hand, was not directly used. The latter enters, however, 
in the definition of the two-point function (^\TX n i(t)X n i(tf)\^f) , where |^) is the coherent ground 
state. To leading order in 1/N and after a Fourier transform one obtains 

Ar.yOO = / dte^-^WTX^Xnjit'm = c% ^ffl^. • (19) 

J S ~r %TJ 

The fact that the residue at the pole has to be 2 NT 2 N (§) = 1, leads to zn = N — 2. 
Besides the Goldstone modes there also exist other excitations. They can be made explicit in di- 
agonalizing the remaining part of This is a straightforward procedure which can be found in 
0. In short, since the non-diagonal part of H 2 is at most bilinear in the operators A\, Af, (3, (3 + , a 
generalized Bogoliubov rotation of the type 

Qt = X u (3 + - Y v p + U U A+ - V U A 1 , (20) 

can be performed and leads to uncoupled modes at the minimum of the action. The diagonalization 
is done by recalling the usual Rowe equations of motion || 

(RPA\ [5Q U , [H 2 , Qt]] \RPA) = n v (RPA\ [5Q U , Q+] \RPA) , (21) 

where \RPA), the full ground state of the theory at this order, is a random-phase approximation 
(RPA) ground state, defined by Q U \RPA) = 0. The Hamiltonian can then be written in the RPA 
phonon basis, \v) = Q^\RPA) , as follows 

N 

H = NH + Erpa + £ n J2 a t a i + E VvQtQv + 0(N- 1 2) . (22) 

i=l i/=±l,±2 

and contains three terms of order (V^/V) 2 , (V^/V) 1 , (\^N)°, respectively. The coefficient of the term 
vanishes. The contribution Erpa = {RPA\H2\RPA) is the RPA correction to the ground-state energy 
and will not be given explicitly here. The frequencies £l u are solutions of the characteristic equation 
of the RPA eigenvalues problem and given by 



r\ 8gs 2 



V = j + 1 ■ ( 23 ) 

In the exact symmetry limit (77 = 0), there exist a pair of zero-energy solutions among the four RPA 

eigenvalues which correspond to two uncorrelated Goldstone modes []. This point is not the mean 
1 Here again, we disregard the infrared problem since the model is only used for demonstration purposes. The reader 
is referred to [tl| for a thorough study of these questions in four space-time dimensions. 



purpose to the present note and therefore will not be discussed further. The reader may consider 
looking into ref.JTJ for a complete treatment of this question. 

We therefore see that the Hamiltonian in Eq. (^) is the same as in ||], however, augmented by the 
'single-pion' term J2iLi a t a i- This extra term arises necessarily in our approach where single bosons 
and pairs of bosons are treated on the same footing. In the single boson state has been treated on 
a heuristic level by neglecting exchange contributions to the self energy. So implicitly, this amounts 
to the same as using ( |22| ) at the order considered. The present systematic scheme puts the treatment 
of ref. |l]] on a firm theoretical ground. 

4 Conclusion 

In this paper we have extended previous work on the Holstein-Primakoff boson expansion for boson 
pairs applied to a relativistic field theory of interacting bosons |]l|] . The aim was to treat simultaneously 
single bosons and pairs of bosons which is necessary to unambiguously define the two-point function 
for the Goldstone mode and to extend the formalism to finite temperature. 

The mapping was applied to the anharmonic oscillator with broken 0(N + l)-symmetry. It was 
explicitly shown that the extension to accommodate single bosons indeed renders, to leading order 
of the 1/A-expansion, N uncoupled Goldstone- as well as RPA phonon modes. This result is novel 
and inaccessible to the bosonic Holstein-Primakoff mapping for boson pairs. The latter is only able to 
provide RPA phonon modes as previously shown in ref. [jlj. The full power of the formalism will reveal 
itself in working out the next-to-leading order of the 1/iV-expansion by providing an unambiguous 
computation of all n-point functions. It also allows for a natural and straightforward extension to 
finite temperature. These two points will be discussed in a forthcoming publication. 
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